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Abstract. One of the remaining issues in the problems of dwarf novae is whether or not enhanced mass transfer 
due to irradiation of the secondary stars could occur during outbursts. In a previous paper (Osaki and Meyer 
2003), we presented a theoretical analysis that shows no appreciable enhancement of the mass outflow rate. This 
conclusion is challenged by Smak (2004) who claims that equations used in our analysis were incorrect and that in 
systems with short orbital periods substantial enhancement could occur. In this letter, we examine the origin of 
such divergent conclusions. We show that Smak's solutions are unacceptable from the standpoint of the equation 
of continuity and that our analysis is an appropriate one to treat this problem. 
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1. Introduction 

It is well known in theory of outburst mechanisms of dwarf 
novae that two rival models (i.e., the disk instability model 
and the mass-transfer burst model) have been competing 
with each other. The disk instability model is now favored 
because of substantial evidence both in observations and 
theory and is widely accepted as the correct mechanism 

' (see, e.g., monographs by Warner 1995, and Hellier 2001). 
However, claims for evidence of enhanced mass overflow 
due to irradiation of the secondary stars during outbursts 

. of dwarf novae, in particular during the superoutbursts 
of SU UMa stars, still appear in the literature from time 
to time. We have critically examined the observational 
evidence (Osaki and Meyer 2003) and concluded that it 
is not well substantiated. Furthermore, we have presented 
a theoretical analysis which shows that irradiation during 
outburst should not affect the mass transfer rate. However, 
this theoretical analysis has recently been challenged by 
Smak (2004) who claims that our equations were incorrect. 
The purpose of this letter is to examine critically these two 
papers and to clarify what the correct approach to treat 
this problem is. 
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2. The problem 



During an outburst of a dwarf nova the surface of the 
secondary star is irradiated by radiation from the cen- 
tral white dwarf, the inner part of the accretion disk, and 
the boundary layer, and the temperature of irradiated re- 
gions of the secondary star is raised. However, the inner 
Lagrangian point from which mass overflow occurs and 
the equatorial region of the secondary star are shielded by 
the accretion disk and only the region of higher latitude 
is exposed. Unequal irradiation then causes fluid flow over 
the surface of the secondary. The question is whether this 
flow is capable to transport heat to the inner Lagrangian 
point and thereby significantly enhance the mass overflow 
rate. 

Controversy has arisen between Osaki and Meyer 
(2003) and Smak (2004) over how to calculate this fluid 
flow. Osaki and Meyer (2003) treated it as a steady 
geostrophic flow using the Eulerian equation of fluid mo- 
tion while Smak (2004) instead calculates trajectories 
of fluid particles for a given initial condition using the 
Lagrangian equation of motion. Although there are several 
minor differences between the two papers such as treat- 
ment of viscosity, coordinate systems used, and temper- 
ature of the irradiated part of the surface, they are not 
essential and thus we concentrate here on the fundamen- 
tal difference in the treatment of the hydrodynamic flow. 
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3. Fluid dynamics in a rotating frame of reference 

The basic equations of fluid flow in a frame of reference 
rotating with the binary system (see Pedlovsky 1982) are 
the equation of motion, 



dv VP T 

— = V* - 20 x v+ — , 

at p p 

and the equation of continuity, 



dp 



+ pV ■ v = 0, 



(1) 



(2) 



where v is the flow velocity in the corotating frame, P 
and p are pressure and density, \I/ is the total gravitational 
("Roche") potential with contributions from both binary 
stars and the centrifugal force, O is the angular velocity of 
rotation of the binary system, T is the viscous force, and 
symbols with bold-face signify vectors. The Lagrangian 
time derivative 4? follows a particular fluid element and is 
related to the Eulerian derivative M- by 



d_ 

It 



d_ 

dt 



+ (v • V). 



(3) 



The equation of continuity is rewritten in Eulerian form 



as 



(pv) = 0. 



(4) 



In unperturbed state in which irradiation is not yet 
applied, we assume that the secondary star rotates syn- 
chronously and is in hydrostatic equilibrium so that 



v =0 



and 



VPp 
Po 



(5) 



where subscript denotes unperturbed equilibrium values. 
4. The two approaches 

We now consider the effect of irradiation. As discussed 
above, the equatorial zone of the secondary star is shad- 
owed by the accretion disk and only a region of higher 
latitude on the secondary star is heated by irradiation. 
The resulting pressure difference between these two re- 
gions drives a fluid flow. Since the depth of the zone af- 
fected by irradiation (which we denote by D) is limited 
to a sub-surface zone a few scale heights deep, we have 
D <C R, where R is the stellar radius. We may then con- 
sider flow in this surface zone only. The flow is essentially 
horizontal, the ratio of the horizontal to the vertical com- 
ponent of velocity is of the order of D/R as seen from the 
continuity equation. The treatment of this flow is com- 
pletely different in Osaki and Meyer (2003) and in Smak 
(2004). We examine each of these two approaches. 



4.1. Osaki and Meyer (2003) 

In Osaki and Meyer (2003) we assume a steady geostrophic 
flow over the surface of the secondary star. Since the pres- 
sure gradient force caused by uneven heating is directed 



in latitudinal direction toward the equator and since in 
geostrophic approximation the Coriolis force balances the 
pressure gradient, the flow is zonal, i.e., directed parallel 
to the equator and, for an inviscid fluid, the equator-ward 
velocity component is zero. Osaki and Meyer (2003) also 
considered effects of a possible turbulent viscosity but con- 
cluded that the equator-ward flow velocity is so slow that 
no effective enhancement of the mass overflow rate is ex- 
pected. 

Osaki and Meyer discussed the particular case of the 
2001 outburst of WZ Sge. Basic quantities there are the 
80 min orbital period of the binary and a half thickness of 
10 9 - 4 cm for the equatorial shadowed belt. When an out- 
burst occurs and the higher latitude region of the sec- 
ondary star is heated by irradiation, a flow sets in. Its 
initial phase may be rather complicated but it will settle 
in an essentially steady state because the typical time- 
scale of an outburst is of the order of a few days, much 
longer than the rotation period of the secondary star. 

The order of magnitude of the left hand side of equa- 
tion (1) is given by 



dv 
~dt 



dv 
~dt 



+ (v • V)v = 



v y2 

t' L 



(6) 



where V, L, and r are characteristic values of velocity, 
length, and time-scale of the flow. The ratio of the inertial 
term to the Coriolis force on the right hand side is the 
so-called Rossby number, Ro, 



rfv/dt 
20 x v 







1 



V 



(7) 



where Ct± is the component of the angular velocity per- 
pendicular to the local surface. 

In our case the Rossby number given by Ro = 
1/(2Q±t) is much less than one except for the equatorial 
region. We may estimate Ro ~ 0.01 for the particular case 
considered by Osaki and Meyer in which Cl± ~ 10 _315 s _1 
and r ~ 10 5 s because we are interested in the flow system 
on the surface on a time-scale of days. In such a case we 
may neglect the inertial term on the left hand side of the 
equation of motion with respect to the Coriolis term. The 
equation of motion for an inviscid fluid then is 



20 x v 



VP 

P 



V*. 



(8) 



Let us use a local Cartesian coordinate system with 
coordinates x in the direction of constant latitude (east- 
ward), y in the direction of constant longitude (north- 
ward), and z normal to the equipotential surface (verti- 
cally upward). Since the main gradients are perpendicular 
to the shadow band, we neglect gradients in the x direc- 
tion and obtain (Osaki and Meyer 2003) 



dP 
dy 



2pv x Qj_ = 0, 
2pv y tt± = 0. 



(9) 
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Without viscosity we have the"geostrophic" flow along the 
shadow band 



1 1 dP 
2f2j_ p dy 



1 Ac 2 
20 1 I 







(10) 



where c = (P/p) 1 ^ 2 is the isothermal sound velocity and 
/ is the half width of the shadow band. Thus the flow 
velocity is at right angle to the driving pressure gradient 
which is balanced by the Coriolis force, similar to flows in 
weather systems and oceans on the rotating earth. Thus 
without friction no heat can be advected into the shadow 
region below Li. The order of magnitude of the zonal flow 
generated by the geostrophic balance in equation (10) is 
V ~ 10 4 cm/s where we used Ac 2 ~ 10 cm 2 /s 2 and L ~ 
10 9 4 cm, which yields a corresponding Rossby number 
Bo = V/(2fl±L) ~ 0.003, confirming the self-consistency 
of the geostrophic flow. 

Here we may add the following notes. As in our pre- 
vious paper (Osaki and Meyer 2003), we used above a 
local Cartesian coordinate system to simplify the com- 
plex Roche geometry. However, the detailed geometry is 
not important in the geostrophic approximation because 
the geostrophic flow is directed perpendicular to the pres- 
sure gradient and as the latter is directed from the shad- 
owed to the irradiated region the flow direction is along 
the shadow boundary. Flow over the surface of a mass- 
losing secondary star in a semi-detached binary was al- 
ready discussed by Lubow and Shu (1975) who argued 
that the horizontal flow is parallel to the isobars calling the 
approximation "astrostrophic" rather than "geostrophic" . 
Our equation (10) is exactly the same as their equation 
(66) in their section V. Oka et al (2002) made numerical 
simulations of the surface flow of a gas-losing secondary 
star and their numerical simulations basically confirmed 
Lubow and Shu's astrostrophic wind. 

4.2. Smak (2004) 

Smak (2004) calculated the motion of fluid elements us- 
ing the Lagrangian equation of motion. He begins the 
calculation at the irradiated side of the shadow bound- 
ary (which we call "the starting point") and then follows 
the motion of each fluid element along its trajectory as 
a one-dimensional initial-value problem. The fundamen- 
tal objection to this approach is that it does not describe 
a continuous hydrodynamic flow but only calculates the 
motion of a free particle that abruptly starts at a par- 
ticular location following the force of a pressure gradient 
obtained from cooling of the flow. Although unphysical, 
let us examine Smak's approach in more detail. 



4.2.1. Quasi-steady case 

Smak (2004) apparently considers a quasi-steady state 
when in his equation (14 ) he calculates the temperature 
gradient from its Lagrangian time derivative: The com- 



plete Lagrangian derivative of temperature T is 
dT _ dT 

dt dt dt ds 

dT dT dT 

dt x " dx* v * dy* ' 



dT dT 
v • VT = — +v- 



(11) 



where s is the distance from the starting point measured 
along the trajectory and x* and y„ are local Cartesian co- 
ordinates used by Smak with cc* directed toward the equa- 
tor and y* directed in longitudinal direction. In Smak's 
formulation the first and third term of the right hand side 
are neglected so that 



dT 
~dt 



dT 

dx* 



(12) 



i.e., steady state is assumed. Accordingly, in his Figure 
2 the trajectories of fluid elements shown are considered 
as stream lines of the flow. The third term on the right 
hand side in equation (|11|) was neglected by Smak because 
he assumed dT/dy* = as his stream lines were fairly 
parallel. 

We therefore adopt the assumption of steady state to 
further examine Smak's solutions. Taking the scalar prod- 
uct of the equation of motion (Eq. (17) of Smak) with the 
velocity v, we obtain 



dw 
~dt 



VP 



V*. 



(13) 



Using the equation of state P — (JStpT/p), with 5ft and 
p gas constant and mean molecular weight, respectively, 
equation l|13|l becomes 



dv 5ftT v ■ Vp 



dt 



H p 



5ft 

—v • VT 

M 



V*. 



(14) 



Smak then assumed that the pressure gradient is 
mainly produced by the temperature gradient so that a 
density gradient is neglected and also takes v • V\I/ = 
because the flow is confined to the equipotential surface. 
Then we obtain 



ldv 2 
2~dT 



5ft dT 
p dt 



(15) 



In Smak's formulation, the Lagrangian time derivative of 
temperature on the right hand side of equation (|15J) is 
given by radiative cooling of the surface element (his equa- 
tion (16)) and he apparently numerically integrated his 
equations as an initial value problem. However, equation 
(fH>jl can be integrated analytically to yield, along a given 
stream line, 



1 



5ftP 
p 



const. 



(16) 



This corresponds to Bernoulli's equation. 

As initial condition Smak apparently assumed that the 
material at the starting point is at rest, v — at s = 0. 
Otherwise one should not have chosen the starting point 
at the shadow boundary. 
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Equation (|16fl then becomes 



2 V 



A* 



KT 

3 

A* 



(17) 



where To is the temperature at the starting point. We note 
that this relation holds for all the trajectories for a given 
set of model parameters listed in Smak's Table 2. 

At first sight the above solution might look reasonable 
because gas at rest at the shadow boundary is accelerated 
to higher velocity by the temperature gradient produced 
by cooling of the flow. However, this solution has the fatal 
flaw that it does not satisfy the equation of continuity: If 
one were to estimate the mass flux from the irradiated to 
the shadowed side by this formula, he would come to the 
unbelievable conclusion of no mass flux across the shadow 
boundary because the velocity there is zero. 

One can demonstrate that Smak's solution of equation 
(16) is applicable only as a steady-state solution for su- 
personic flow. However, as in Smak's solution matter near 
the shadow boundary is subsonic, it can not be reconciled 
with the equation of continuity (Eq. 0J which in steady 
state is 



v • V/9 + pV ■ v = 0. 



(18) 



A subsonic flow behaves more like an incompressible fluid 
and the first term on the left hand side of equation (|18J) 
becomes negligible, v • Vp = 0, an approximation used by 
Smak in expressing pressure gradient by temperature gra- 
dient in his equation (13). However, Smak's accelerating 
flow has a non-zero velocity divergence in contradiction to 
the equation of continuity fEa. I18JI . In particular, in his 
solution the velocity divergence becomes infinite at the 
shadow boundary from where gas at rest is accelerated to 
a finite velocity: In Smak's form, equation (|15fl is written 
as 

dv x . VtdT 



jji dt 



In the shadow zone the right hand side of equation (|19|l is 
non-zero so that dv Xt jdx* becomes infinite when v Xr = 0. 

This decisive defect of Smak's solution results from 
his assumption that matter at the shadow boundary stays 
at rest. But there is no steady flow on parallel equator- 
ward steamlines with zero velocity at the boundary. If 
instead gas has a finite velocity across the shadow bound- 
ary, matter from higher latitude must move in to replace 
matter that moved into the shadow region. Since this mat- 
ter has already a significant zonal velocity when passing 
the shadow boundary, a picture very different from that 
of Smak emerges. 

One may ask whether at the shadow boundary cool 
material might well up from below to fill the gap created 
by gas drained equatorwards, and always replace it by 
matter with zero zonal velocity. However, in order for up- 
welling gas to prevent inflow from higher latitude it must 
provide the same pressure at the same depth as the neigh- 
boring irradiated region. With vertical hydrostatic equilib- 
rium this means the same temperature down to the same 



depth. To heat up a column by irradiation to that depth 
requires the same time as that for which the existing hot 
material was exposed i.e. the time since the start of irra- 
diation. As soon as this time exceeds the time on which 
heated matter flows in, of the order of the sound travel 
time (about an hour), up- welling of gas is suppressed by 
inflow. 

4.2.2. Unsteady case 

This shows that Smak's approach is unacceptable as a 
description of steady flow. Let us examine whether his 
solutions might reflect an unsteady initial flow. In this case 
we suddenly switch on the pressure gradient and observe 
what happens. As Smak argued, the matter at the heated 
side of the shadow boundary first moves along the pressure 
gradient, i.e. equator- wards. This lowers the density at 
the shadow boundary and creates a pressure deficit with 
respect to neighboring higher latitude regions from which 
gas then moves in to fill the void. Since this gas has lower 
specific angular momentum than that which already left 
the starting point, it has there a significant zonal velocity. 
This velocity increases with time as material from higher 
and higher latitude is drained into the shadow region. This 
adjustment occurs over all the surface of the secondary 
star until the flow at the boundary of the shadow region 
has a zonal component sufficient for the Coriolis force to 
balance the pressure gradient, and our steady geostrophic 
state is reached. 

One may estimate the time needed to establish 
geostrophic equilibrium. In the case of WZ Sge discussed 
above and for the minimal model of a spherical secondary 
star already gas from a distance of 10 9 6 cm from the equa- 
tor, on reaching the interior of the shadow belt would 
have enough zonal velocity to balance pressure gradient 
by Coriolis force. The hydrodynamic time for this matter 
to arrive in the shadow belt is again only of the order of an 
hour, the same order of magnitude as the orbital period. 

5. Conclusion 

We conclude that in dwarf novae outbursts enhanced mass 
transfer by irradiation of the secondary star is prevented 
by the strong Coriolis force. 
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